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Abstract 



We resolve a discrepancy between the 577(2) spacial string tension at finite temper- 
ature, and the value obtained by monopoles in the maximum Abelian gauge. Previous 
work had incorrectly omitted a term due to Dirac sheets. When this term is included, 
the monopole and full SU(2) determinations of the spacial string tension agree to within 
the statistical errors of the monopole calculation. 



The string tension in pure £77(2) lattice gauge theory has recently been calculated 
from monopoles, giving results that agree with the full SU(2) string tension, to within 
the statistical errors of the monopole calculations |J. Given this success, it is disap- 
pointing that at finite temperature, monopoles apparently fail to explain the fundamental 
representation spacial string tension, a s . In the region of temperatures above the decon- 
fining transition, the values for a s obained from monopoles are too small, by an amount 
well outside of statistical errors, and the discrepancy increases with increasing tempera- 
ture 1,1. 

In this paper we resolve this discrepancy. In the part of the calculation in which 
Wilson loops are calculated from monopoles, a term involving Dirac sheet variables had 
been dropped. When this term is included as it should be, the monopole and full SU (2) 
determinations of the fundamental spacial string tensions agree to within statistical errors. 
The present calculations were carried out on lattices of size 16 3 x N t a.t (3 = 2.5115, for 
N t = 4, 6, 8, 12. At this value of /?, the deconfining transition corresponds to N t = 8 [|J, 
so N t = 4, 6 are in the high-temperature phase. We compare to the full SU(2) results for 
a s obtained by Bali et al || at this same (3 value on 32 3 x N t lattices. Before presenting 
our results, we briefly review the way the monopole calculations are done, and discuss the 
correction due to Dirac sheets. 

A summary of the steps |I[ involved in the monopole calculations is as follows: (1) 
SU(2) configurations are projected into the maximum Abelian gauge. (2) The resulting 
SU(2) links are factored into Abelian and charged parts. (3) The Abelian link angles are 
used to locate the magnetic current of the monopoles. (4) Monopole Wilson loops are 
calculated. (5) The string tension is extracted from fits to the monopole Wilson loops. 
Steps (3)-(5) proceed exactly as they would in U{1) lattice gauge theory. 

For SU (2) lattice gauge theory, these steps involve strong assumptions. Briefly, these 
are that in the maximum Abelian gauge (MAG), long range physics is Abelian and con- 
trolled by monopoles. Each link of a configuration in the MAG is approximated by its 
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Abelian part, exp(z0^r 3 ), where r 3 = 03/2 is the isospin generator, and — 2n < <^ < 2ir. 
The key assumption for Wilson loops is 

< W SU (2) > ~ < >, 

where the ~ sign means equivalent at long range. The Abelian Wilson loop is given 
by 

W4 = exp(i 4>n{ x ) J »{ x ))- 00 

In Eq.(||), M is a rescaled link angle, M = 4>l/2, while J M is an integer- valued current 
describing the path traversed by the heavy quark. Motivated by results in £7(1) lattice 
gauge theory ||, the Abelian Wilson loop is further assumed to factor into a short 
range part involving the exchange of neutral gluons, times a term W mon arising from 
monopoles. Thus for the confining part of the heavy quark potential 

< W sui 2) > ~ < W$ > ~ < W mon > (2) 

is supposed to hold. In particular, Wsu(2), W^, and W mon should all produce the full SU(2) 
string tension. Here we are concerned with the correct calculation of W mon , and therefore 
with the second equivalence of Eq.(^j), < > ~ < W mon >. The assumption of Abelian 
dominance and the use of the MAG are mainly concerned with the first equivalence in 
Eq.(|). 

To calculate W mon , monopole variables must be located in a M configuration J7|. 

Abelian plaquettes 4>^ u are formed from the link angles M , and expressed as (p^ u = 

<fi' + 2nm^ u , where <f>' G [— ir, ir], and m^ u is an integer. The surface formed by the 

dual variables m*^ describes the Dirac sheets which are present. [] Two currents can be 
1 Our conventions are that direct lattice link variables (j^, At), originate from the labelling site, while 
dual lattice link variables (m^, A™) terminate at the the labelling site. Direct lattice plaquette variables 
(m M „, F^ v , F™*) have the site at the lower left corner of the plaquette, while dual lattice plaquette 
variables (m*„, F^*, F™) have the site at the upper right corner. In going from the direct to the dual 
lattice, the discrete difference <9^ is replaced by d^, and vice versa. 

3 



derived from m^, one the magnetic current of the monopoles, m M = d^m*, the other 
the electric Dirac sheet current, = d~m^ v . Geometrically, m M flows on the edge of the 
open surface defined by m*^ u . To obtain W mon , we start from an expression which can be 
derived analytically in (Villain) U(l) lattice gauge theory ||, 

W mon = exp(2vrz £ J^Afa)). (3) 

X,fl 

Here is an electric vector potential whose source is j M . The subtlety with which the 
present paper is concerned arises when an attempt is made to express W mon in terms of 
the magnetic current m M . This cannot be done using vector potentials, since the magnetic 
current produces a magnetic vector potential A™, whereas the electric current of the quark 
can only couple to an electric vector potential, in this case A e . Progress can be made 
by going to field strengths. We introduce a surface D^ u of plaquettes whose boundary is 
the heavy quark current J M , so that = d~D^ v . Then using the lattice form of Stoke's 
theorem, W mon can be written as the exponential of a flux integral 

2iri 

W mon = exp(— £ 2V(*W*)), ( 4 ) 

X,fl,U 

where the field strength is given by F^ u = d^Al — d*A e , and satisfies the electric Maxwell 
equation, d~F^ u = j^. Analogous to A e , there is a magnetic vector potential A™, whose 
source is m M . The field strength given by = d~A™ — d~A™, satisfies the magnetic 
Maxwell equation, d^F™ = m^. 

In a finite volume, it is not valid to replace F* in the exponent of Eq. (f|) by the dual 
of F™, defined as usual by F™ = ^e^pF™^. For the case of periodic boundary condition 
in all directions, the correct equation relating F* v and F™* is 

F; u {x) + F™*(x) = m^(x) - m^, (5) 

where rh^ is the space-time average of m^, 

™^ = 77Y, m A x )- ( 6 ) 

* X 
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It is straightforward to derive Eq.([|). Here we simply list various consistency checks. 
Applying d~ to Eq.(|5|) gives an identity since d~F™* = , and d~F^ u = d~m^ u = j^. 
Likewise, the equation obtained by applying d+ to the dual of Eq. (^|) is identically satisfied. 
The constant —fn^ on the right side of Eq.(|5|) is needed to make the space-time average 
of the right hand side vanish. The space-time average of the left hand side vanishes for 
periodic boundary conditions, since F^ u and F™* are linear combinations of gradients in 
their respective vector potentials, so the sum over all x of + F™* vanishes. 

Using Eq.(||) to replace in Eq.@, we have 

27T? 

W mon = exp(-— J2 D, v (x)(F™(x) + m, v )). (7) 

The integer m^ v term on the right hand side of Eq.(|5|) has no effect on W mon . However, 
the presence of the non-integer term in Eq.(^) means that in addition to the magnetic 
current, the six numbers must be specified to obtain W mon . Since Eq.(^) is equivalent 
to Eq.([|), and therefore to Eq.([|), the necessary property that W mon is unity for a plane- 
filling loop is now guaranteed. This is false if the fh^ v term is omitted in the exponent of 
Eq.(|7|). While these three forms for W mon are completely equivalent, there are reasons for 
preferring the representation of Eq.(^). If monopoles are the cause of non-perturbative 
phenomena, it is desirable to express physical quantities in terms of to the greatest 
extent possible. Further, the magnetic current is sparse, occupying only a small percentage 
of the links of the lattice, and is essentially limited to the values ±1. The electric Dirac 
sheet current, j M , on the other hand, is dense, occupying a large fraction of the links of the 
lattice, and is not dominated by the values ±1. Finally, varies under local deformations 
of Dirac sheets. The deformation — > m M „ + d^n u — <9+n M , an integer, leaves m M 
and m M „ unchanged, while j M does change, j M — > + d£(d~ ■ n) — (<9 + • <9~)n M . Although 
W mon is invariant under local Dirac sheet deformations, there remains a dependence on 
the overall Dirac sheet topology. For example, a loop of magnetic current could have a 
Dirac sheet consisting of plaquettes with non- vanishing m*^, which tile the "inner" area 
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<?.(SU(2)) 
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0.049(3) 


0.065(3) 


0.0643(6) 


6 


0.026(1) 


0.040(4) 


0.0381(4) 


8 


0.022(2) 


0.031(3) 


0.0325(7) 


12 


0.029(1) 


0.034(1) 





Table I 



Table 1: Spacial string tensions from magnetic current, magnetic current + Dirac sheets, 
and for full SU(2) 

of the loop, or by virtue of the periodic boundary conditions, the "outer" area. The 
difference between the value of W mon for these two cases comes entirely from m^. 

We now turn to our results. The methods we used for generating SU(2) configurations 
and our gauge- fixing criterion are described in detail in [1]]. After equilibration, we saved 
configurations every 20 lattice updates, resulting in a total of 500 configurations of mag- 
netic current m M , and the six Dirac sheet space-time averages, fh^ ||. As we previously 
found at zero temperature, the magnetic current is sparse. The average fraction of links 
carrying magnetic current can be conveniently written as / x 10~ 2 . For our lattices of 
size 16 3 x N t the values of / were 1.30(1), 1.07(1), 1.18(1), 1.25(1), for N t = 4, 6, 8, 12 re- 
spectively. The averages of Triply £1X6 statistically zero, so as a measure of the effect of 
in Eq.(0), we give 2ir times the standard deviation of rh^, over the three purely spacial 
planes of the lattice. Expressing this as h x 10~ 3 , we have h = 11(2), 8(2) 9(2) , 5(1), for 
N t = 4,6,8,12. 

Monopole Wilson loops were calculated using Eq.(0). These loops are located in the 
purely spacial planes of the 16 3 x N t lattice. One dimension of the loop can be regarded 
as a separation R, while the other S is a pseudo-time, and by the same arguments as used 



6 



in the usual case of a symmetric lattice, W mon (R, S) should approach exp(—SV ps (R)), 
as S becomes large, where V ps (R) is the pseudo-potential. The values of V ps (R) were 
determined by fitting — \n(W mon ) to a straight line in S, over the interval S m i n = R + 2 
to S max = 13. Then linear-plus- Coulomb fits were performed on V ps (R) over the interval 
R = 2 to R = 7. The coefficient of the linear term in these fits gives our estimate of the 
spacial string tension due to monopoles. We denote the spacial string tension deduced 
from W mon values obtained from Eq.(|7]) as cr s (m M , m^). For our old results where the 
m M „ term was omitted in Eq.(|7|) ||], we use the symbol ^(m^). In table I , we compare 
these two determinations of a s to the full SU(2) results of || on 32 3 x N t lattices. As can 
be seen by a glance at the table, the a s (m^, rh^ u ) values agree well with the full SU{2) 
results, whereas the values obtained by omitting the contribution due to clearly 
too small. 

The values of the fh^ are of course highly correlated with the magnetic current in 
a given configuration. However, the fh^ cannot be constructed from knowledge of the 
magnetic current alone. A useful equation relating fh^ to F™* can be derived by summing 
Eq.(H) over the \i — v plane. It is possible to show that fh^ v can be recovered from F™* 
using this equation only if | J2x m fj,u\ < |A*/3, where A a p is the area in lattice units of the 
plane dual to the // — v plane. If this restriction is ignored and the resulting estimate of 
fn^y is used in Eq.(|7|), values of o s which are intermediate between those of columns 1 
and 2 of Table I are obtained. We have also checked that using the term by itself to 
calculate W mon does not produce an area law for Wilson loops, so the effect of rh^ is not 
a simple additive term in the string tension. 

In summary, the apparent serious discrepancy between monopole and full SU(2) an- 
swers for the spacial string tension has been removed. The price paid is that a rather large 
contribution from Dirac sheets must be included. Dirac sheets are no longer "invisible" 
in a finite volume. It will clearly be of great interest to further explore this effect as a 
function of spacial volume and temperature, including zero temperature. 
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We take one MAG gauge-fixed configuration for each original 577(2) configuration. 
Based on recent work by A. Hart and M. Teper, (to be published) our conclusion that 
inclusion of Dirac sheets effects is essential is likely to survive in a more elaborate 
calculation involving several gauge copies for each configuration. 
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